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ABSTRACT 

We use the eight year light-curve database from the MACHO (MAssive Compact Halo Objects) 
project together with infrared colors and magnitudes from 2MASS (the Two Micron All Sky Survey) 
to identify a sample of 22,000 long period variables in the Large Magellanic Cloud (referred to hereafter 
as LMC LPVs). A period-luminosity diagram of these stars reveals six well-defined sequences, in 
substantial agreement with previous analyses of samples from OGLE (Optical Gravitational Lensing 
Experiment). In our analysis we identify analogues to galactic LPVs in the LMC LPV sample. We find 
that carbon-dominated AGB stars populate only two of the sequences, one of which includes the Mira 
variables. The high-luminosity end of the same two sequences are also the location of the only stars 
with J — K s > 2, indicating that they are enshrouded in dust. The unknown mechanism that drives 
the variability of stars in the longest period produces different morphology in the period-luminosity 
diagram as compared to the shortest period sequences, which are thought to be caused by pulsation. 
In particular, the longest period sequence extends to lower luminosity RGB stars and the luminosity 
function does not peak among the AGB stars. We point out several features which will constrain new 
models of the period-luminosity sequences. 

Subject headings: galaxies: individual (LMC) — stars: AGB and post- AGB — stars: variables: other 



1. INTRODUCTION 

Long period variables (LPVs) are giant stars that pul- 
sate with periods ranging from weeks to years. The pro- 
totypical case is Mira (Omicron) Ceti, from which all 
large-amplitude (> 2.5 magnitudes in V) regular LPVs 
take their name. LPVs are typically asymptotic giant 
branch (AGB) stars that are in the late stages of stellar 
evolution, where the pulsation mechanism(s) - and the 
relationships between pulsation, mass loss, and the even- 
tual ejection of the stellar envelope as the star becomes 
a planetary nebula - are all poorly understood. During 
this relatively short AGB phase, stars manage to lose 
a significant fraction of their total mass, returning pro- 
cessed material to the interstellar medium. A qualitative 
picture of the mass-loss process invokes instabilities in 
the helium and hydrogen burning shells surrounding the 
degenerate core of the star. This provides the energy for 
a variety of stellar pulsation modes, causing shocks in the 
envelope which lift material high in the atmosphere. The 
temperature is then low enough to allow the formation of 
dust, which is subsequently ejected from the atmosphere 
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by radiation pressure from the luminous core. Pulsa- 
tions - which affect both the radius and temperature of 
the star - together with the shocks and dust formation, 
combine to produce significant photometric variability. 

The study of LPVs has undergone a recent revival 
with the advent of surveys producing large catalogs 
of variable stars. Gravitational m i cro-len sing surveys, 
such as OGLE - iPaczvnski et all lll994fl and OGLE 
II - lUdalski. Kubiak. fc Szvmanskil (Il997f>: EROS - 
lAubourg et aU 11199,51): M ACHO - IMcock et al.l Ipcjl : 
and MOA - iBond et all l|200lD require that millions of 
stars be imaged nightly. The variable star catalogs that 
arise as a by-product contain thousands of stars with 
light curves co vering severa l thou sand days. 

As noted bv iCook et al.1 l)1996j) the long period vari- 
ables in the MACHO database for m five para l lel se- 
quences in period-luminosity space. iWood et aU l)1999f) 
classified these and suggested possible underlying mech- 
anisms for the five. Our use of infrared magnitudes in 
the period-luminosity relations shown in Figure ^ re- 
veals six sequences. Wood's second sequence ("B" in 
his notation) is split into our Seq uences 2 and 3. In 
etgrms of our observed sequences, IWood et alJ l|1999fl 
found that Sequence 1 corresponds to Mira type vari- 
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ables, which are theorized to be pulsating in the fun- 
damental mode. Higher-order pulsations are invoked to 
explain the shorter period sequences (2, 3, and 4). The 
sequence that lies at longer-periods and lower luminosi- 
ties than the Miras, E, showed light curves that sug- 
gested that all of these variable stars were eclipsing bi- 
naries. Stars in the longest period sequence, Sequence D, 
proved more mysterious since the long secondary periods 
of these stars could not be explained by radial pulsations 
- they are longer than the periods of Miras, which them- 
selves are theorized to be pulsating in the fundamental 
mode. Additionally, they all exhibited multiple periods 
with the short period associated with Wood's Sequence 
"B". This prompted Wood et al. to propose that Se- 
quence D was caused by the eclipse of the AGB star by 
a cloud of material around an unseen compani on. 

More recent studies usin g OG LE -II data lit a et al.1 
(12001: iKiss k Bedding! ll200l; llta et all (|200" l 7 
Kiss k Beddinel (I200J7 : ISoszvnski et alJ l)2004j) : 
Groeneweg cn (2004J)] have extended the work on 
the period-luminosity sequences. These authors showed 
that the use of infrared magnitudes (e.g. from 2MASS) 
for luminosity split the original second sequence into 
two sequences. Further, their data show that Wood's 
first two sequences (our first three sequences) exhibit a 
split in luminosity due to the contribution of red giant 
branch (RGB) stars. Most previous work on LMC LPVs 
overplot each perio d of these predominate ly multiperi- 
odic stars. However. IKiss fc Bedding (|2004fl find there is 
no noticeable change in the period-luminosity sequences 
when plotting only the strongest period, or even when 
plotting the second strongest. For this reason results 
can be compared between those works and o nes which 
find on ly a single peri od for each star: e.g. llta et all 
( 2002). llta et al] l|2004h and this paper. 

In this paper we use the full eight-year MACHO 
database to re-examine the LPV period-luminosity re- 
lations. MACHO observations have several advantages 
over previous studies: the time baseline of the light 
curves is twice as long (eight years compared to four 
years), and the colors of MACHO light curves yield su- 
perior sampling of photometric variability of RGB stars. 
As with previous studies, we use 2MASS infrared data 
to supplement our optical lightcurves. 

2. DATA 

2.1. The MACHO LMC Variable Star Database 

The MACHO project ijAlcock et al.l 11997ft comprises 
eight years of observations of the Large and Small Mag- 
ellanic Clouds and the Milky Way Bulge. In this work 
we se lect LMC stars from the MACHO variable star cat- 
alog ijAlcock et al.l l2003'). Stars from the full database of 
several million stars were selected for this catalog if the 
central 80 percent of points in the light curve failed to fit 
a constant magnitude in a %-squared test. This criterion 
resulted in 207,632 variable candidates in the LMC. 1 

MACHO emp loyed a nonparametric phasing technique 
l)Riemannlll994D to find periods for each candidate vari- 
able star. Data were taken simultaneously in a red and 
a blue filter, and were analyzed independently to find a 
period and amplitude for each variable candidate. There 

1 The MACHO variable star catalog is online at: 
http : //wwwmacho . monaster . ca/ 



is no substantial difference between the periods found, 
but a plot of the red versus the blue amplitudes shows a 
bias toward larger amplitudes in the blue. We adopt the 
blue periods and amplitudes so that the wider spread 
in amplitude will allow us to better discriminate types 
of LMC LPVs. The period-finding algorithm tends to 
alias noisy data and stars with chaotic variability. Thus, 
stars with period aliases at the total survey length, one 
year, and multiples of one day, in particular up to the 
fourth multiple and down to one-ninth (|) of a day, were 
all removed from the catalog. Note that these alias cuts 
result in blank vertical stripes in the period-luminosity 
diagram (Figure Only 52 percent of the LPVs have 
well-determined periods, a total of 21,441 stars. We are 
currently analyzing the ~ 20, 000 LPVs with poorly- 
determined periods with the goal of understanding the 
multi-periodic nature of these stars. 

The amplitudes in the MACHO database are deter- 
mined from the difference between the median of five 
points nearest the maximum and five points nearest the 
minimum in the light curve. A model light curve is used 
to identify the times of maximum and minimum light, 
so stars that are not well fit by the model will tend to 
average to zero amplitude. Thus, noisy or chaotic data 
will not be assigned meaningful amplitudes. 

2.2. 2MASS Photometry 

The Two Micron All Sky Survey (2MASS) measured 
J, H, and K s magnitudes for a pproximately half a bil- 
lion objects over the entire sky l|Cutri et al.N 2003). The 
sky coverage, resolution, and sensitivity of 2MASS make 
it a very useful source of infrared magnitudes. We have 
matched the MACHO and 2MASS data by position for 
LMC LPVs with well-determined periods. Since the 
2MASS magnitudes were taken at random phases we ex- 
pect scatter in the K magnitudes from the intrinsic vari- 
ability of these LPVs. Although most of the stars in our 
sample are small amplitude pulsators and contribute a 
negligible amount of scatter in the K s band, some of the 
stars that make up Sequence 4 are Mira variables. Mi- 
ras are the largest am plitude LPVs and have a typical 
K amplitude of - 0.6 (IWoodl I2000D . Thus the 2MASS 
magnitudes of the largest amplitude stars may deviate 
from the mean by as much as 0.3 mags. We accept this 
as a source of scatter, and, assuming all LMC stars are 
at the same distance, we use the 2MASS K s magnitude 
as our primary luminosity indicator. 

Infrared magnitudes can also be used to distinguish 
the state of late AGB stars. After evolving through the 
RGB and horizontal branch, stars begin to ascend the 
AGB with a high O/C ratio in their atmospheres. These 
oxygen-rich AGB stars include both the so-called early 
AGB stars and AGB stars that have had their first ther- 
mal pulses. Low mass stars will remain on the AGB long 
enough that thermal pulses initiate the "third dredge- 
up" , the deepening of their convection zones down into 
material enriched in 12 C. The spectral type then changes 
to a carbon star. The obj ective-prism survey of the LMC 
bv iKontizas et all l)2001l) showed that carbon stars ap- 
pear on the infrared color-color diagram in a red tail 
running from the stellar locus. Infrared color cuts of 
J - H > 0.89 and H - K s > 0.32 isolated 83% of 
the carbon stars in their sample, with only a nine per- 
cent contamination rate. A single color cut in J — K s 
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Fig. 1. — LMC period-luminosity diagram showing the stars with well-determined periods. The luminosity split at the top of the RGB 
(K s = 12.3 ± 0.1) is observed in sequences 2, 3, and 4. 



can also select these stars since it is perpendicular to 
the tail of carbon stars; the result is a sample that is 
less complete but also has less contamination. In addi- 
tion, the 2MASS color-magn itude diagram of the LMC 
( Nikol aev fc Weinberg! l200fy shows a clear division at 
J — K s = 1.4 due to carbon stars. In this paper we will 
draw particular attention to LPVs near this color bound- 
ary with the caveat that we are missing some twenty per- 
cent of the carbon stars - those that are warm and hence 
b luer than J — K K = 1.4. 

iNikolaev fc Weinberg! l)2000|) also find a population of 
very red stars (J — K s > 2) that lie along the reddening 
vector from the carbon stars. They propose that these 
stars are surrounded by a dusty circumstellar envelope. 
Both AGB stars and at least a few protostars are known 
to lie in this region. Although such dusty LPV are likely 
to be very late AGB stars, pin-pointing the exact evo- 
lutionary state of these latter stars will require a better 
understanding of the connection between mass-loss, pul- 
sation, and the late stages of AGB star evolution. 

3. DISCUSSION 

Period-luminosity relations for variable stars in the 
LMC are shown in Figure The two Cepheid se- 
quences are clearly visible between < log P < 1 and 
13 < K s < 16. When plotted in K s , as shown here, 
it is evident that Wood's original sequences split in both 
period and luminosity. It also appears that the low lumi- 
nosity ends of the three shortest period sequences are off- 
set to longer periods. These effects have been previously 
seen in OGLE II data: lKiss fc Bedding! l(2003l ). lita et all 
(2004). The lower luminosity stars in Sequences 2, 3, 
and 4 have been identified as RGB stars ijKiss fc Bedding! 
120041) using an analysis of the second derivative of the 
luminosity function. The number of stars located above 
and below the tip of the RGB in each pulsation sequence 
is given in Tabic ^ For the tip of the RGB w e use 
K s = 12.3 ± 0.1 from INikolaev fc Weinberg! l|2000t) who 
analyzed the luminosity function of all 2MASS stars in 
the direction of the LMC. Note that some previous au- 
thors have used K s = 12.1, found by analyzing the lumi- 
nosity function of only the variable stars. 

3.1. Connections Between LMC and Galactic LPVs 

Galactic LPVs are identified primarily in The Gen- 
eral Catalog of Variable Stars ijKholopov et al.l Il996j) 
which has traditionally classified LPVs by their V band 
behavior. The classification scheme, tabulated in Ta- 
ble |2 separates LPVs into Mira types and classes of 
semi-regular (SR) variables. Mira variables are de- 
fined as those pulsating with a period greater than 80 
days and an amplitude of at least 2.5 in V. Whereas 
Mira variables show regular and strong periodicity, one 
would imagine that "semi-regular" variables would have 
weak or poorly defined periods. Indeed, SRb stars have 
poorly expressed periodicity, or multiple periods, or 



TABLE 1 

Properties of the Period-Luminosity Sequences 



Sequence Star Count Amplitude 





above TRGB a 


below TRGB a 


Mean b 


Max 


4 


1337 


673 


0.11 


1.75 


3 


2780 


540 


0.17 


1.90 


2 


3967 


310 


0.34 


3.79 


1 


5009 


628 


0.77 


8.04 


E 


636 


1384 


0.28 


4.47 


D 


2163 


1802 


0.38 


5.58 



a Tip of the Red Giant Br anch; K s = 12.3 ± 0.1 in the LMC 
I Nikolaev fc Weinber J 2000 ). Above the tip of the RGB the LPVs 
are all AGB stars. 

b We are unable to estimate a minimum amplitude of pulsation as 
MACHO's amplitude findi ng al gorithm reports zero amplitude for 
noisy or chaotic data (see Ij2.ll . 



only occasional periodicity, or chaotic pulsation. How- 
ever, SRa stars are not necessarily semi-regular at all; 
they are similar to Miras but with sm aller amplitudes 
ijLebzelter. Schultheis. fc Melchio'rll2002|) . 

Due to the empirical nature of the GCVS classification 
system, connections drawn to the observed LMC LPV 
period-luminosity sequences are often n ot simple or di- 
rect. H owever, connections can be made: iLebzelter et alJ 
(2002) find that Miras in the Galaxy and in the Mag- 
ellanic Clouds have the same period-luminosity rela- 
tionship, implying that the pulsation mechanism is 
not strongly dependent on mctallicity. Also, parallel 
period-luminosity sequenc es exist in Galactic samples 
ijBedding fc Ziilstral 1)19981) : iKnapp et all <|2003j) L as in 
the LMC. 

Figure shows LMC LPVs of different amplitudes 
located in period-luminosity space, and Table ^ gives 
the mean and maximum pulsation amplitudes for each 
sequence. There is only a small difference between 
our Bmacho amplitudes and the V band amplitudes 
used t o classify Galactic variables. Bes sell fc Germany! 
( 1999) found transformations from MACHO magnitudes 
to Cousins V and R. For B M acho — Rmacho > 1-0 

V — Bmacho — —0.07 — 0.10(Bmacho — Rmacho)- 

The change in the amplitude is thus due to color 
changes during the star's variability cycle. In particular, 
the change in the amplitude between V and Bmacho 
is restricted to 10% of the color change in the MACHO 
photometric system. 2 

With reference to the LMC LPVs in Figure Se- 
quence 1 is composed on the long period side by Mi- 

2 Since we kn ow the typical color of AGB stars, Bmacho ~ 
Rmacho ~ 1-5 lAlcock et al. 1999), we note that V magnitudes 
will be ~ 0.22 brighter than the Bmacho magnitudes. 
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Fig. 2. — Period-Luminosity diagrams highlighting stars in seven amplitude bins (with divisions at Bmacho = 0.07, 0.17, 0.3, 0.7, 
1, and 2.5). The background distribution (in gray) is for all of the LMC LPVs from FigureQand is shown as a reference. Note that the 
highest amplitude stars are the Mira variables, located on the long period side of Sequence 1. 



Fig. 3. — Period-Luminosity diagrams highlighting stars in five J — K a color bins (with divisions at J — K s = 1.1, 1.2, 1.4, and 2). The 
background distribution (in gray) is for all of the LMC LPVs from Figureffland is shown as a reference. Note that stars with J — K B > 1.4 
are carbon stars, and stars with J — K a > 2 are obscured by dust (see H2.2I . The bluest stars (J — K a < 1.1) show the RGB populations 
of the period-luminosity sequences. 



ras. iCioni et al.l <)2001[) found that higher amplitude SRa 
stars lie on Sequence 1 with the Miras; low amplitude 
SRas, on the other hand, lie on the shorter period se- 
quences (Sequences 2, 3, and 4). ICioni et all (|2001f) in- 
vestigate stars from Sequence D and find that they are all 
multiply-periodic SRb stars. iWood et alJ l)1999l) found 
that the shorter period of these stars usually fell on his 
sequence "B" , which corresponds to our Sequences 2 and 
3. Wood's work did not identify any stars from Sequence 
D with a primary period on Sequence 4. Unlike sequence 
1, the highest amplitude stars fall on the shorter period 
side of the sequence. 

The RGB stars at the low luminosity end of Sequences 
3 and 4 are low amplitude pulsators, and a similar set 
of stars is visible at the base of the Sequence D. The 
latter sequence's stars are particularly interesting as their 
variability mechanism remains unknown. Sequence D is 
also notable for a significant number of high amplitude 
variables (1 < B amp < 2.5). 

3.2. Color and Luminosity Functions 

Figure shows LMC LPVs of different J - K s color 
located on the period-luminosity diagram; as discussed in 
H2.2I stars with J — K s > 1.4 are likely to be carbon stars, 
while stars with J — K s > 2 are likely obscured by dust. 
Interestingly, the most heavily obscured stars lie only at 
the high luminosity end of Sequence 1 and are Miras. The 
majority of carbon stars lie at the high luminosity end 
of Sequences 1 and 2. Sequence D shows that the lower 
luminosity, small amplitude pulsators visible in Figure [2] 
also have bluer colors than the rest of the sequence, again 
indicating a separate RGB population. 

The color-magnitude diagrams and luminosity func- 
tions for each sequence are shown in Figure the lu- 
minosity functions are scaled such that the relative pop- 
ulation of stars in each sequence can be compared. The 
background distribution in the color-magnitude diagram 
is for the entire MACHO variable star catalog, whereas 
the background lumino sity function is for all LMC stars 
measured by 2MASS l|Nikolaev fc Weinberdl2000(h 

The "bumps" in the luminosity function for each se- 
quence correspond to particular populations at differ- 
ent stages of stellar evolution. The top of the RGB is 
the largest such bump, the tip of which was found by 
INikola,ev fc Weinberg! pWil) to be K s = 12.3 ± 0.1. Se- 
quences 3 and 4 show a substantial contribution from the 
top of the RGB, whereas Sequences 1 and 2 are composed 
primarily of stars on the AGB (which overlaps with the 
RGB at this magnitude). Sequences 1, 2, 3, and 4 all 
show a peak in the luminosity function at K s ~ 11.3; 
this peak in variable AGB stars is not seen in the overall 



LMC luminosity function ijNikolaev fe Weinberd feOOO) . 

The types of stars in each sequence are obtained by 
comparison with the color-magn i tude diagram of the 
LMC from iNikolaev fc Weinberg] IpOOO]) . Sequences 3 
and 4 have similar distributions including a population 
at the tip of the RGB as well as oxygen-rich AGB stars. 
The stars at J — K s » 1.25 with K s < 10.25 are young, 
massive AGB stars; they are most commonly found in 
Sequence 3. Sequences 1 and 2 also have similar pop- 
ulations, with both oxygen- rich AGB stars and a large 
contribution from carbon stars. Sequence 1 extends fur- 
ther to the red in J — K s , probably due to dust ob- 
scuration around those stars. Sequence E is dominated 
by RGB stars, as can be seen from both its luminos- 
ity function and position in the color-magnitude dia- 
gram. This is expected if it is composed of eclipsing 
binaries. Both Sequence D and E lack high luminos- 
ity young, massive AGB stars and the AGB dominated 
population that produces the significant peak in the lu- 
minosity function at K s ~ 11.3 in Sequences 1, 2, 3, and 
4. Interestingly, Sequence D ends abruptly at K s ~ 13.7 
(Sequence E is artificially cut off at K s = 14.5). There 
is no feature in the "deep" LMC color-magnitude dia- 
gram (Ni kolaev fc Weinberg! 1200(1) associated with this 
position, or with the bottom of the first four sequences. 
The difference in the low-luminosity cut-off of Sequence 
D compared with the first four sequences, as well as the 
lack of a bump in the luminosity function at K s ~ 11.3 
in Sequence D, suggests that the variability is caused by 
a different mechanism than the radial pulsation modes 
proposed to explain Sequences 1, 2, 3, and 4. 

3.3. Comparisons to Theory 

Explanations of the structure in the period-luminosity 
diagram often invoke radial pulsations, with parallel se- 
quences indicating populations of overtone pulsators. 
We have plotted the linear, rad i al, no n-adiabatic pul- 
sation models of IWood fc Sebol l|199(iD with our LMC 
LPV data in Figure These models were forced to 
fit the observed LM C Mira period-luminosity relation of 
iFeast et all (^89). Note that the observed Mira rela- 
tion lies precisely where we see the largest amplitude 
pulsators in Figure [3 Sequences 1, 2, 3, and 4 have 
all been previously attributed to overtone pulsations, 
but Sequence 4 seems to lie apart from the rest. Se- 
quences 3 and 4 share every important feature discussed 
in this paper - morphology in the color-magnitude di- 
agram, luminosity function, infrared colors, and, to a 
lesser degree, amplitudes (see Table H). New models are 
required in order to decide whether Sequence 4 repre- 
sents the next highest overtone pulsation, or if there are 
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TABLE 2 

GCVS Long Period Variable Classification 



Type 


Spectral Type 




Period (days) 


Amplitude (V) 


P-L sequence (s) 


Mira 


late a giants 


80 - 


1000 


> 2.5 


4 


SRa 


late a giants 


35 - 


1200 


< 2.5 


1,2,3,4 


SRb 


late a giants 


20 - 


2300 




1,2,3,6 


SRc b 


late supergiants 


30 - 


"several thousand" 


~ 1 




SRd b 


F, G, & K 


30 - 


1100 


0.1 - 4 





a a "late" giant specifically indicates a spectral type of M, C, S, Me, Ce, or Se. 
b These rare types are not represented in our collection of late giants. 



Fig. 4. — Color-magnitude diagrams and luminosity functions (along the left vertical axes) for each period-luminosity sequence arranged 
from shortest to longest period (see Figure 0. The background distribution in the color-magnitude diagra m is for the entire MACHO 
variable star catalog, whereas the background luminosity function is for all LMC stars measured by 2MASS (Nikolacv & Weinberg 2000). 
The tip of the RGB is marked at K s = 12.3 in each panel. 



Fig. 5. — Period-luminosity diagra m over laid with the ob- 
served LMC Mira relation of Feast ct al. (1989) (dashed line) and 
the 3rd, 2nd, and 1s t overtone (respectively, from left) models of 
IWood Kr, Sebol I1996N (solid lines). The fun damental mode of the 
models was forced to fit the Mira relation of Feast ct al. 1 1989). 



damped pulsation modes between Sequences 3 and 4. Se- 
q uence 4 may not be shor test period-luminosity sequence 
as lSoszvnski et aD l)2004D find evidence of a sparcely pop- 
ulated sequence at still shorter period. 

Sequ ence D was originally theorized by IWood et alJ 
( 1999) to be caused by eclipses from a dust-enshrouded 
unseen companion. This hypothesis required that ~ 
25% of AGB stars exist in se mi-detached binaries. 
IWood. Olivier, fc Kawalerl l)2004D have recently com- 
pleted an analysis of four years of echelle data for three 
Galactic stars with lon g secondary periods (examples of 
which are tabulated in Houk (1963)). Additional photo- 
metric data for 111 stars from MACHO allow them to in- 
vestigate, and rule out, explanations involving solely ra- 
dial pulsations, non-radial pulsations, orbiting compan- 
ions, non-spherically symmetric stars, dust obscuration, 
and chromospheric activity. They propose that the most 
likely explanation is a low-degree g + mode in an extraor- 
dinarily thick radiative layer, which would allow large 
amplitudes at the stell ar surface, comb i ned w ith large- 
scale star spot activity. Sos zvnski et alJ l)2004|) find that 
Sequence D is made up of two populations: stars that 
do not have one of their periods on sequence 4 and stars 
that do. This split seperates Miras and SR variables from 
very small amplitude red variables respectively. Thus Se- 
quence D is composed of stars which cannot be uniquely 
assigned to any specific class (AGB, RGB, Mira, SRa, 
SRb, very small amplitude pulsators, etc.). 

4. CONCLUSIONS AND FUTURE WORK 

The combination of observations over a long time base- 
line, accurate optical photometry extending below the tip 
of the RGB, and 2MASS infrared photometry has allowed 
us to place multiple constraints on the stars inhabiting 
each of the sequences observed in the period-luminosity 



diagram of the LMC. 

• Most of the LPVs in this analysis, which comprises 
the 52% of LMC LPVs with well-determined pe- 
riods, are SRa stars or Miras. The amplitude cri- 
terion (SV > 2.5 for Miras) used to differentiate 
Galactic SRa and Mira stars seems arbitrary, as 
SRa's and Miras overlap in our diagrams. Stars 
from Sequence D are SRb stars as previous work 
has found that they are never monoperiodic. 

• Our period-luminosity Sequences 3 and 4 both 
show low luminosity extensions comprised of RGB 
stars, and have AGB stars that are oxygen rich 
(J — K s < 1.4). As discussed in t )2.2l this implies 
that these AGB stars haven't yet undergone the 
third dredge-up. Compared with the significant 
population of redder (carbon) stars in Sequences 
1 and 2 this implies that older stars are segregated 
to the latter two sequences. The populations of 
young, massive AGB stars in Sequences 3 and 4 
that are identified in Figure 0] (see 13.2(1 support 
this hypothesis. 

• The shortest period sequence (Sequence 4) is signif- 
icantly seperated in period from the longer period 
sequences but it is otherwise undistinguished. 

• Sequences D and E tend to follow the LMC lumi- 
nosity function more closely than the others, indi- 
cating that they are more nearly drawn from the 
general population of giant stars in the LMC. In 
particular, neither sequence shows a "bump" in its 
luminosity function due to AGB stars, unlike Se- 
quences 1, 2, 3, and 4. These characteristics sug- 
gest that a different mechanism produces the vari- 
ability of these sequences. Sequence E is likely pop- 
ulated by eclipsing binary systems. 

• Although the physical mechanism that causes Se- 
quence D is still elusive, t here are numerous con- 
straints on possible models. ISoszvnski et al.l |2004) 
show that this sequence is composed of bright high 
amplitude Mira and SR variable stars and a dim- 
mer population of small amplitude red variables. 
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This latter population is bluer, extends to lower 
luminosity, and is broader in period than the RGB 
star components of Sequences 2, 3, and 4.t The 
AGB population in Sequence D exhibits higher am- 
plitude pulsation on its shorter period side, oppo- 
site what is found for Sequence 1. 

We are currently analyzing the ~ 20, 000 multiply pe- 
riodic LMC LPVs that are not included in the present 
sample, with the goal of understanding the relationships 
between the period-luminosity sequences and the posi- 
tions in period-luminosity space of stars at various stages 
of their evolution. 
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